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We present an integral transformation capable of extracting the moments of arbitrary Paley- 
Wiener entire functions f{z) against a given spectral distribution dP{u>) from the short-time values 
of the correlation function C{t). The transformation reads 



f{u)dP{uj) 



1 i^mCH^)^^ 



ds. 



It is proven to be valid for correlation functions that can be extended analytically to the entire 
complex plane with the possible exception of two branch cuts on the imaginary axis. If the analytic 
continuation exists only for a strip |Im(z;)| < ro, then the integral transformation remains valid 
for all Paley- Wiener functions obtained by Fourier-Laplace transforming a compactly supported 
distribution, with the support included in the interval (— 2ro,2To). Finally, if the support of the 
distribution is contained in the interval (— To,ro), then the generalized moment can be evaluated 
from the short-time values of the correlation function exponentially fast. 
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I. INTRODUCTION 



Several authors have investigated the possibility of reconstructing spectral 
functions from imaginary-time correlation functions by inverting a real Laplace 
transfor m^i^''^i^i^i^i^i^i^ or solving a Hamburger moment problem.^'^"'^^ It is assumed that 
the input information is represented either by the correlation function C (t) for imaginary 
values of i in a limited interval or, more strongly, by the values of C{t) in some open disk 
about the origin in the complex plane. That the information in imaginary-time is greatly 
supplemented by its real-time counterpart, in as much as complementary aspects of the 
spectrum are described, has been recognized a while ago by Kim, Doll, and Freeman^ as 
well as Krilov and Berne»2, My own effort to produce a path integral technique capable 
of furnishing information in an entire open diskii has been motivated not only by the 
physical insight obtained from the aforementioned research, but also by an intriguing 
observation of Lynessi^ that some problems that are numerically unstable on the real 
line become numerically stable when formulated in the complex plane. Lyness' most 
explicit example was the problem of numerical differentiation, which he proved to be 
stable if performed by means of Cauchy's integral formula. 

Although the specific autocorrelation functions vary, they can always be casted in the 
form of a Fourier transform of a probability measure (see, for example, Ref. ■IL'l 

C{t) = [ e''^'dP{uj). (1) 

In some formulations, the actual measure can be signed, but the negative and positive 
parts are always finite and, by linearity, no loss of generality is incurred by assuming 
dP{ijj) to be a positive measure. C{t) is also called the characteristic function of dP{uj), 
whereas C{—it), if it exists, constitutes the moment generating function. In many cases 
of practical importance, the imaginary-time correlation function C{—it) is finite for t in 
some range (— ro,To). tq is usually proportional to h,(3 and is rarely infinite. Because of 
the absolute integrability, it can easily be argued that C{t) is in fact analytic in the entire 
strip with |Im(z)| < tq. Baym and Mermin^^ were perhaps the first to point out that 
this well-known fact in statistics automatically implies that the correlation function can 
be reconstructed from its imaginary-time values by analytic continuation. The latter 
quantities are normally assumed to be easier to obtain, for example by path integral 
Monte Carlo techniques. 
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Analytic continuation is, of course, a numerically unstable procedure. The quality 
of the reconstructed spectral functions depends tremendously on the quality and the 
"form" in which the input information is given. In this paper, we assume that the 
information is available in the form of the values of the correlation function in an open 
disk {z € C : |z| < tq} about the origin, a disk included in the region of analyticity of 
C{t) (that is, only short-time information is available). By virtue of Lyness' observation, 
a Hamburger moment problem can be set up, with the moments of the spectral function 
computed by means of a contour integral. If 7 is a Jordan curve about the origin and 
situated in the region of analyticity of C{t), we obtain 



4t) 



t=o 



fc! 
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C{-iz) 



dz. 



(2) 



Taglianii^ has employed this same observation to good effect in order to substitute an 
inverse Hamburger moment problem for an inverse real Laplace transform. Given the 
known instability of the moment problem when monomials are utilized a natural 
question that arises is whether or not we can compute moments of more suitable func- 
tions, for example polynomials that are orthogonal with respect to some other measure 
or more general complex entire functions. It is the purpose of the present work to show 
that this is possible in several cases. 



II. COMPUTING GENERALIZED MOMENTS OF SPECTRAL FUNCTIONS 

I have been experimenting for some time with a modification of Eq. (jJJ of the form 

Eq. (jSJ is equivalent to Eq. (O because fc! — s^e'^'ds. Nevertheless, it allows for im- 
proved accuracy in evaluating moments of polynomials that are orthogonal with respect 
to another probability distribution, perhaps one that closely resembles dP{uj). Indeed, 
if f{z) = X)fe=o o-kZ^ : then one can define the function 

and use linearity to justify the equality 

e-'Vf{s)ds^ [ f{uj)dP{uj). (5) 
Jr 

The catch here is that the polynomial f{z) may actually be known in the factorized form, 
which is numerically stable. Then, the Gauss-Laguerre quadrature rule of appropriate 
order can be utilized to compute moments of arbitrary polynomials in a stable fashion. 
Notice that Vf{s) is a polynomial of degree n as well and is exactly integrated by the 
Gauss-Laguerre quadrature rule if the number of knots is greater than [n/2] + 1. Modified 
moments against orthogonal polynomials represent a significantly less redundant way 
of storing information. As pointed out by Sack and Donovauj-I they may enable the 
computation of the quadrature knots and weights of the Gauss quadrature with weight 
function dP{u)). This quadrature rule defines a measure made up of point masses that 
can be shown to converge weakly to dP{uo) (see, e.g., Ref. and Theorem 3.12 in 
Chapter 2 of Ref.|18j). 

Nevertheless, it is tempting to utilize Eq. together with more general analytic func- 
tions that are entire on C. For example, a symmetric positive function f[z) that decreases 
to the right and left on the real line, such as the Gaussian f{z) = exp(— z^/2)/-\/27r, con- 
stitutes the perfect example for constructing delta sequences that might provide the 



ds. 



(3) 
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value of the spectrum at a point of continuity. Starting with such a function, we can 
obtain a complete description of the spectrum in the form 

Uf( ^) dn-i (6) 



for any level of resolution cr > desired. Surely, we would like to utilize entire functions 
f{z) with tails that decay fast on the real line. In the following section, we will show 
that this is always possible if the measure dP(uj) is compactly supported. 

Unfortunately, for the more interesting cases, the measure is not compactly supported. 
Its tail only decays as an exponential of the form exp(— tqIljI), for some tq > 0. A 
suggestive example is given by the Miller, Schwartz, and Tromp flux-flux spectral function 
for the free particle^^ which has the continuous distribution 

where Ki{x) denotes the respective modified Bessel function of the second kind. Setting 
To = hp/2, the tail of the distribution behaves roughly as jwl"""^/^ exp(— To|a;|). The 
consequence is that the distribution does not integrate exponentials of order greater 
than or equal to tq and, therefore, the correlation function is not defined for purely 
imaginary times of the form {ir : \t\ > tq}. The actual correlation function reads 

CO = J-^^lZHll-, (8) 

and is seen to exhibit two branch cuts on the imaginary axis, branch cuts that form the 
set {z G C : z = it, \t\ > tq}. 

Nevertheless, it is the case with this correlation function and, perhaps, with many other 
correlation functions of practical importance that they are analytic everywhere with the 
possible exception of two branch cuts. Under this assumption, in Section IV, we prove 
that 'Pf{s) is integrable against e~^ds and that Eq. ||SJ| is true for all entire functions 
of Paley- Wiener type. These functions are given by the Fourier-Laplace transforms of 
some compactly supported distribution /(k), i.e., 

f{z) = I f{K)e'-^d^. (9) 

In the most general case, the correlation function can be extended only to a strip of 
analyticity defined by the inequality |Im(z)| < tq. In such cases, Eq. © is valid for those 
Paley- Wiener entire functions which are the Fourier-Laplace transform of distributions 
with the support included in the compact interval (— 2to, 2to). 

Having established the physical and mathematical motivation, we dedicate the remain- 
der of the paper to the actual proofs. 



III. THE COMPACTLY SUPPORTED CASE 



Theorem 1 LetdP{uj) be a probability measure onM such that its characteristic function 
is analytic at the origin and let tq > be the radius of convergence of the Taylor series 
about the origin (tq may be infinite). Let f : C C be an entire analytic function, the 
real restriction of which increases slower than any exponential, i.e., for any e > there 
exists a constant M > such that < Me*^'^', V x G M. Then Vf{s) is an entire 

function. If dP(uj) is compactly supported, then 

e-'Vf{s)ds= [ f{cj)dP{u;). (10) 
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Proof. By hypothesis, the Taylor series of / about the origin f{z) = X]fc>o '^fe-^'^ 
is absolutely convergent on the whole complex plane, whereas the Taylor series of the 
moment generating function C{—it) — X]/£>o Mfe^*^/^' convergent in an open disk of 
radius tq. It follows that the contour integral expressed by Eq. Q is convergent for all 
s S C. By Cauchy's theorem and the uniform convergence of the partial sums on the 
contour 7, the result is the convolution series 



k=0 



(11) 



'Pf{s) is an entire function since it is a power series convergent for all s € C. 

We treat the remainder of the theorem by starting with some preparatory material. 
The measure dP{uj) integrates to a finite number any function of exponential order 
r < To- Indeed, 

\f{cj)\dP{Lj) <M [ e^l"ldP(w) <M [ (e™ + e-™)dP(w) 

= M[C{iT) + Ci-ir)] < 00. (12) 

Thus, the right-hand side of Eq. H1U|) is finite. 
Let (7 : C ^ C be defined by 



fc=0 



A:!' 



(13) 



g{z) is also entire and, for any e > 0, there is > such that \g{z)\ < Meexp(e|z| 
These assertions are established at once by the sequence of inequalities 
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(14) 



^k>0 



and the finitude of the last series in parenthesis (the series, which is a power series, is finite 
by the absolute convergence of the Taylor series of / about the origin for an arbitrary 
radius). In fact, Eq. H14|) also establishes that the absolute series J2k>o Wk\\z\''/kl is of 
exponential order e, for all e > 0. Thus, the series g{ujs) = J2k>o'^k{i^s)'' /kl can be 
integrated term by term against both measures e~'^ds and dP{u>). We obtain 



.fc=0 



dP{uj) 



g{ujs)dP{uj) (15) 



and 



00 00 1 " /" 

/(^) = E"'='^'' = E"'='^''i:t / 
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ds 



s^e-'ds 



e ^g{u)s)ds. 



(16) 



The reader may already see why Eq. (|10|l could be true. It boils down to proving a 
change in the order of integration. Indeed, by Eq. (|16|l . 



f{uo)dP{u) = 



e ^g{ujs)ds 



dP(w), 



(17) 
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whereas by Eq. H15|) and assuming Vf is integrable against e "ds, 



/•oo /"OO r /• 

/ e-"7'/(s)ds= / e-" / 



g{ujs)dP{uj) 



ds. (18) 



The change in the order of integration as weU as the integrabiUty of 7^/(s) could be 
justified by the Fubini-Tonelli theorem provided that the kernel g{i-us) were integrable 
against the product measure. This is true for measures dP(uj) supported in some arbi- 
trary interval [— r, r], with r > 0. Since in this case \^k\ < 2r'^, we calculate 

/•oo /• oo oo 

/ e-'ds / dP{Lo)\g[ujs)\ < ^ \ak\\fik\ < 2^ |afe|r'= < ex., (19) 

•^0 k=0 k=Q 

again by the absolute convergence of the Taylor series of f(z). The proof of the theorem 
is completed. □ 
A quick look at Eq. (|19l) shows that the ease with which we have demonstrated the 
validity of Eq. H1Q(I in such a general setting is due to the convergence of the series 

oo oo oo 

F{e) = J2we'' = Y,ak6'' / uj'^dPiu:) = / Y,^k{Oujf dP{uj) 

k=0 k=0 k=Q 

f{euj)dP{uj), yoec. (20) 



That the Taylor scries of f{9z) can be integrated term by term follows from the fact 
that the domain of integration is some finite interval [— r, r]. Thus, the right-hand side 
of Eq. (|20|l is always finite because f{9Lo) is bounded on [— r, r], since it is continuous. 
Then, the uniform convergence of the Taylor series justifies the change of order between 
integration and summation and establishes the finitude of the left-hand side series. 

If the measure dP{uj) is not compactly supported, we cannot expect the left-hand 
side series to be convergent for all 9 & C, unless the function /(z) is rather trivial 
(e.g., a polynomial). In fact, the series may not be convergent at all for the more 
interesting functions, such as f{z) — exp(— z^) = J2k>o(~^'^)'^ Indeed, if dP{uj) = 
exp(— |aj|)/2(iti), the even moments of which are /i2fe ~ (2fc)!, we compute 

F{9) = Y.-,^,^o'-^li~^)'^^o"'^ 

k=0 k=0 

which is divergent for all 7^ (by the ratio test). On the other hand, if f{z) = exp(iz) = 
S/c>o(*^)'^/^'' interesting thing happens. The series 

00 oo> 

F(0) = 5]a,Mfce'==^(*e)2'= = — ^, 1^1 <1, (21) 

k=0 k=0 

is convergent if \9\ < 1, but divergent if 1611 > 1. Nevertheless, the integral 

Fi9) = jj{9u)dP{u,) = i^e-l-le^^-dc. = |lm(0)| < 1, (22) 

is convergent for all 9 of imaginary part in the interval (—1, 1). 

Clearly, in this case, F{9) is just the correlation function of the measure dP{uj). In 
the introduction, we have already mentioned that F(9) is analytic in the strip defined 
by the equation |Im(6')| < tq (for the example described by Eq. tq = 1). Given the 
relevance of the convergence of the series expansion of F{9) in the proof of Theorem ^ 
it makes sense to restrict our attention to entire functions that exhibit such properties. 
With respect to their ability of recovering information via delta sequences, our intuition 
says that functions like f{z) — exp(— z^) are, perhaps, of little help. Functions like 
f{z) = sin(z)/z, which exhibit properties similar to f{z) = exp(iz), may actually do the 
job. In the next section, we show that the latter assertion is true, indeed, under the 
assumption that the correlation function of dP{uj) can be analytically continued to the 
whole complex plane save two branch cuts on the imaginary axis. 
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IV. THE GENERAL CASE 



For the remainder of the present paper, we shaU assume that the function f{z) and 
the measure dP{uj) satisfy the conditions of Theorem ^ save the compactness of the 
support of the measure. In fact, the functions f{z) will be restricted to a smaller subset, 
in agreement with the opinion expressed in the last paragraph of the preceding section. 
The class is provided by the following definition. 

Definition 1 An entire function f : C ^ C is said to be of Paley-Wiener type (or 
simply Paley- Wiener function ) if it is the Fourier- Laplace transform of some compactly 
supported distribution /(k), i.e., 

f{z) = I f{K)e'-^dK. (23) 



The distributional version of the Paley-Wiener theorem asserts that if the distribution 
is supported in the interval [— with i? > 0, there exist positive constants C and 
N such that 

|/(z)|<C(l + |z|)^e^|i-WI. (24) 

(The converse is also true). Thus, the function f{z) has at most a polynomial growth for 
real arguments. In fact, a Paley-Wiener function f{z) has at most a polynomial growth 
in any strip |Im(2:)| < r, with t > 0. 

Proposition 1 Let f : C C be a Paley- Wiener function and let S be the strip S = 
{9 eC: \lm{9)\ < to/B}. Then, 

[ \f{euj)\dP{u;) <oo, \/9eS. (25) 
Jr 

The function F{9) : — > C defined by 

F{e) = I f{euj)dP{u;) (26) 
Jm 

is analytic. More specifically, the Taylor series of f{9uj) (regarded as a function of 9) 
about any real number 9r can be integrated term by term and produces the Taylor series 
of F{9) about 9^, a series having a radius of convergence of at least tq/B. 

Proof. We have already mentioned during the proof of Theorem ^ that the measure 
dP(ijj) integrates all functions of exponential order less than tq (according to Ea. I12|l. 
Thus, according to Eq. 1)24(1 . the integral in Eq. (|25|l is absolutely convergent for all 9 
with i?|lm(0)| < To. Even more, if if is a compact set included in the strip S, we have 
J3|Im(0)| < Ti < To and \9\ < T2 for some ti,T2 > and for all 9 £ K [by continuity, 
lm(0) and \9\ attain their extrema at some points in K, a compact set which is properly 
included in S] . Then the inequality 



\f{9uj)\dP{uj) <C / (l + |0w|)^e^'i™(^)'"dP(cj) 

<C J {1 + T2|w|)^e^i"(iP(w) < oo (27) 

tells us that the integral defined by Eq. 1251) is not only finite, but also uniformly bounded 
on each compact set. By the continuity of the integrand and the dominated convergence 
theorem, it follows that F{9) is continuous on S. 

Let 6*^ € R be arbitrary and let f''^\z) denote the /c-th derivative of f{z). Since f{z) 
is an entire function, the Taylor series of f{9u) regarded as a function of 9 

oo ^ 

/(M = ^/W(0,,c.)i^(f?-0,)'= (28) 

fc=0 
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is convergent on the entire complex plane, for all €E K. By Cauchy's integral formula 
(let 7 be a Jordan curve encircling the point 9r), 



f 



dz. 



so that the Taylor series now reads 

oo 



k=0 



2ni 



dz. 



(29) 



Given an arbitrary 9 such that \9 — 9r\ < tq/ B (with 9^ G M, as in the hypothesis of 
the problem), we choose the contour 7 to be a circle centered about 9r and of radius 
T > 0, with r such that |0 — 0r| < t < tq/B. Then the integrability term by term against 
the measure dP{u!) of the series given by Eq. H29|l follows from the Fubini theorem once 
we prove the absolute integrability against the product measure (of dP{uj), the counting 
measure, and the Lebesgue measure for the polar parameterization of 7). Using Tonelli's 
theorem to justify the arbitrary order of integration, we compute 



00 „ 



d(t) 



\9-9r\'' \f[{9r+Te''^)uo\\ 



27r 



oi0 



|fe+l 



1 

27rT 



^fc=o 
1 



2tT {T^\9^9r\) Jo 



( f dP{ou)\f[{9r + Te'^)uj 
dP{oj) \f [{9r + Te">')uj] 



The geometric series is convergent because \9 — 9r\ < t. On the other hand, since 
|lm(0r + re**^)] < r < tq/B, the integral against dP{uj) is convergent for all (p € [0,27r). 
In fact, given that the circle constituting the contour is a compact set included in S, 
the integral against dP{uj), as a function of 4>, is uniformly bounded. It follows that 
the integral against dcj) is finite and the term-by-term integration of the Taylor series is 
established. We therefore have 



00 r „ 

F{9)^Y1 f^''\drUj)uj''dP{uj) 

J n L-^M. 



(9 - 9rf 
kl 



(30) 



Since 9r GM. and 9 with \9 — 9r\ < tq/B are otherwise arbitrary, the proof is concluded. 
□ 

The preceding proposition tells us that the Paley- Wiener functions borrow many of 
the properties of f(z) — exp(iz). The latter function generates the correlation function 
itself 



C{9) 



e'''^dP[uj). 



Thus, we obtain once again that the correlation function is analytic on the entire strip 
|Im(6')| < To (in this case B = 1). The function F{9) can be expressed in terms of the 
correlation function as follows: 



m 



f{K)e"''^'^dK 



dP{uj) 



'dP{Lu) 



dK= J{n)C{9n)dK. 



(31) 

If \9\ < Tq/B, the formal change in the order of integration can be rationalized by means 
of Eq. l|5n)l . as applied for exp{iz) and /(z) about the origin. Indeed, the series 



00 r 

C{9k) = Vk^ / fuj^dP{uj) 
j=0 L^M 



9^ 
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converges absolutely for all \k9\ < tq. As a function of k, the series converges in the 
C°° topology on the open interval |k| < tq/\9\ (the case with 6* = is trivially true). 
This open interval contains the support of the distribution /(k), because \9\ < tq/B. By 
continuity of the distribution as a linear functional, 



/ f{K)C{9K)dK^Y^ [ f(K)K^dK [ i^UJ^dP{, 



9^ 



E 



/W(0)M!dP(.)' 



The order of the summation and integration can be changed again by means of Eq. H30|l . 
as applied to f{z) about the origin. Then the interior series is seen to sum to f{9Lu) and 
its integral equals F{9), indeed. 

We can extend the validity of Eq. I|31() from the disk 16*1 < tq/B to the strip S by 
analytic continuation. Indeed, the left-hand side has been proven to be analytic in the 
strip S. It will coincide with the right-hand side expression on S provided that the latter 
is analytic in 9 on S. The right-hand side expression is well-defined because C(9k), as 
a function of k, is analytic in the open set (regarded as a subset in R) defined by the 
equation |Ini(0K)| = |K||Im(6')| < tq. If 9 ^ S, then this open set contains the support 
of /(k) because |lm(0)| < tq/B. By a similar argument, the first order derivative of the 
right-hand side expression, namely f{K)KC'^^\9K)dK, also exists for all 9^3. 

To continue, let us notice that it is enough to show that, for each 9 ^ S, there exists a 
compact set K Cz S that contains all the complex points of the form zk, with z in some 
neighborhood U of 9 and k in some precompact open set containing [—B,B]. In this 
case, the second derivative of C{zk) against z, i.e., the function C^^^^ {zk)k^ , is bounded 
as absolute value by some constant Af > (by continuity). If 6'„ 9, the Taylor series 
with remainder says that, for n large enough. 



C{9k) - C{9nK) 



C^^'>{9k)k 



- 9 ™P 



C^^^zk)k^ \9 



< 



M. 



Thus, the convergence [C{9k) — C(0„k)]/(6' — 6'„) C'^^\9k)k is uniform with respect 
to K in an open set containing the support of /(k). However, a similar argument with 
C{9k) replaced by C^^\9k)9^ shows that, for every j > 1, 



d^ C{9k) 



dni 



C(6'„k) d^ 
9„ ^ dnJ 



C^^\9k)k 



uniformly on the same precompact open set containing the support of f{n). It follows 
that the convergence [C{9K) — C{9nn)]/ {9—9n) C'^^^{9k)k happens in the C°° topology 
on that set and therefore 



lim 



~ C{9k)-C{9^ . 



C{9k) - C{9nK) 



dn 



= / f(n)C'^^\en)KdK. 



Since the last quantity exists, we obtain that Jg^f(K)C{9K)dK is differentiable at any 
point 9 (z S. It therefore defines an analytic function on S, which must coincide with 
F{9). 

Finding such a compact set K for a given e S* is not difficult. Let e = [tq/B — 
|Im(6l)|]/2. Of course, e > 0, since 9 & S. Because < B{tq/B - e) < B{tq/B - e/2), 
there is 77 > such that {B + ri){To/B — e) < B{to/B — e/2). Consider the set K made 
up of the numbers of the form nz with n G {—B — ri,B + ri) and — 0| < e. Any element 
Kz of K satisfies the inequalities |Im(Kz)| — |K||Im(z)| < {B + ri){To/B — e) < rg — Be/2 
as well as |Re(K2:)| = |K||Re(2;)| < {B + 77)[|Re(6')| -I- e]. If we replace K with its closure, 
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the bounds transform from strict inequalities to inequalities, but the set is still properly 
contained in S. By the discussion in the preceding paragraph, Eq. H31|l is valid on S. 

Nevertheless, the same technique can be applied to analytically continue F{9) on a 
larger connected domain that includes the open disk \0\ < tq/B. Indeed, there was 
nothing special about S except for the existence of an appropriate compact set K asso- 
ciated with a point 9. Assume now that C{d) can be extended analytically on the whole 
complex plane except perhaps for the subset {z G C : z = iy, \y\ > tq}. Then for any 
e eC with rj = Re{9) ^ 0, the set K = {zn € C : k, e {~B ~ l,B + 1), \z ~ 9\ < |?7/2|} 
is included in the domain of analyticity of C (9) . Indeed, one notices that the inequality 
|Re(zK)| = |K|Re(2:) > |k|?7/2 prohibits the set K from intersecting the forbidden region. 
Since K is trivially bounded, its closure is compact. By the continuity of the function 
Re(z), the preceding inequality remains valid on the closure of K, perhaps with the sign 
> replaced by >. It follows that F{9), as now re-defined by Eq. ((HlJ, is analytic on 
the open set Ke{9) ^ 0. We summarize the findings of the last two paragraphs in the 
following proposition. 

Proposition 2 On S = {z € C : |Im(z)| < tq/B}. the following equality of analytic 
functions holds 

F{9) = [ f{n)C{9n)dK. (32) 

In addition, if C{9) can he analytically continued to the whole complex plane except for 
the set {z ^ 'C : z — iy, \y\ > Tq}, then F{9) can be analytically continued to the whole 
complex plane except for the set {z <C : z = iy , \y\ > tq/ B}, by means of Eq. yS^) . 

The relation between 'Pf{s) and F{9) can easily be established from the Taylor series 
expansion of Vf^s), which is given by Eq. If 7 is any Jordan curve encircling the 

origin and enclosed in the disk about origin of radius tq/B, then the Cauchy integral 
formula and the uniform convergence of the Taylor series show that 




T \fe=0 



(33) 

However, after the summation of the series is done, the contour can be extended to the 
strip S or to the whole domain of analyticity of F. 

We now have all the tools to demonstrate the first theorem of the present section. 

Theorem 2 Assume f is a Paley- Wiener function and dPicu) is a probability measure, 
the correlation function of which can be analytically continued to the whole complex plane 
except for the set {z G C : z = iy,\y\ > tq}. Then the restriction of'Pf{s) to [0,oo) is 
integrable against e~^ds and 







'Viis)ds = / f{oj)dPiLo). (34) 
Jm 

Proof. From the definition of F{9), it follows the Eq. (|34|) is equivalent to 



/•oo 

/ e-^Vf{s)ds^F{l) 
Jo 



or, by Eq. to 



e 







1 

2ni 



z~^F(z)e''^dz 



1 



ds = F{1). (35) 



We can take the contour 7 to be the circle centered about 1/2 and of radius r = [1 + 
[tq/ BYy-/'^ /2. Indeed, the circle contains the origin in its interior because its radius is 
strictly greater than 1/2. It is included in the domain of analyticity of F{9) because the 
radius is strictly smaller than [1/4+ (ro/i?)^]^/^, the minimal distance to the forbidden 
set. 
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If z = x + 1/2 + ij/ e 7, we have a;^ + = and {x + 1/2)^ + = + x + 1/4 < 
+ r + 1/4 = (r + 1/2)^. As such, 

zy (x + 1/2)2 + y2 r2 + a; + l/4 - (2r + l)2 ^ 

If M > is an upper abound for |2~^F(z)| on the circle 7, then 

-s{l-l/z)-l. 



F{z)\ < A/cxp{-,sRc(l - z-i)} < Af cxp 



is uniformly bounded by a function that is integrable against the product measure. By 
Fubini-Tonelli's theorem, it follows that F{s) is integrable against e~''ds. By inverting 
the order of integration, we compute 



J— / z^^Flz) 



The proof of the theorem is concluded. □ 
An important corollary of Theorem [21 is the following statement asserting that the 
correlation function can be reconstructed pointwise. In this case, the functions f{z) 
have the form exp(izt) and are obtained by Fourier transforming point masses (Dirac 
functions) about n = t. Setting P/(s) = 'P[f]{s) for clarity and using * to denote the 
dummy variable z, we have 

Corollary 1 In the conditions of Theorem\^ 

/>oo 

C{t) = / e~'V [e"*] {s)ds, (36) 
Jo 

for all t e M. 

Looking back at the proof of Theorem |21 we can see that the reason we require the 
existence of an analytic continuation of F{9) is that the quantity Re (l — becomes 
negative in a circle about 1/2 of radius 1/2. Thus, the Fubini-Tonelli theorem cannot be 
utilized if the contour 7 passes through the disk delimited by this circle. However, if the 
circle of radius 1/2 centered about 1/2 is properly contained in S, and this happens if 
1/2 < tq/B, we can always pick a circle centered about 1/2 of radius strictly greater than 
1/2 that is still properly included in S. This circle contains the origin in its interior. On 
it, the function je"^*^^"^/^^ | is bounded by exp(— es) for some e > 0, whereas |z~-'^_F(z)| is 
bounded by some M > 0. Because the circle is contained in S, the analytic continuation 
of F{9) outside S is no longer needed. However, because B cannot be greater than 
2to, wc lose the ability of constructing delta sequences by rcscaling the argument of the 
function f{z) (the information is band-limited). We obtain the following result. 

Theorem 3 Regardless of whether an analytic continuation of the correlation function 
outside the strip S exists or not, Eq. holds for all Paley- Wiener functions f with 

the property that the support of f is contained in the interval (—2tq,2to). 



V. A LOOK AT THE NUMERICAL IMPLEMENTATION 

In this section, we show that if the support of the distribution /(k) is further restricted 
to be contained in the interval (— tq, tq), then the generalized moment f{uj)dP(Lu) can 
be evaluated exponentially fast from the short-time correlation function. To motivate 
the further restriction of the support, let us start by putting together Eqs. Q and (O 
in a single formula 

/ f{u;)dP{u;)= He-' 



Zm F \z/ z 



ds. 



(37) 
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We have explicitly utilized the square parenthesis in order to emphasize that the right- 
hand side must be regarded as an iterated integral, if the results from the preceding 
sections are to hold. The contour integral in the square parenthesis, that is the value of 
'Pf{s), is calculated by quadrating an highly oscillatory integrand, a scenario that might 
lead to numerical difficulties. For fixed s, results of Fornarc^*^ show that the quantity 
P/(s) can be approximated exponentially fast by numerical quadrature. If we pick the 
contour 7 to be a circle of radius r < tq centered about the origin, we obtain 



Jo 



-ire 



)d(j). (38) 



According to Fornaro's Corollary 5.1 (the corollary will be stated below), the one- 
dimensional integral in Eq. (|38|l can be approximated by trapezoidal quadrature ex- 
ponentially fast. Nevertheless, the rate of convergence is different for different values of 
s and generally worsens in the limit s — > cx). 

Not surprisingly, the negative scenario exposed in the preceding paragraph can be 
avoided provided that the Paley- Wiener functions are further restricted in such a way 
that the iterated integral appearing in Eq. 1)3711 is also valid as a double integral (i.e., 
the conditions of the Fubini-Tonelli theorem are satisfied). In fact, in slightly more 
restrictive conditions, we show that Fornaro's corollary implies that the left-hand side 
of Eq. H37|) can be approximated exponentially fast from the short-time values of the 
correlation function. Let us assume that f{z) is a Paley- Wiener function obtained from 
a distribution /(k) with support included in the interval {—tq,tq). Thus, B < tq and 
we take r such that B < t < tq. From the Paley- Wiener theorem, we obtain that 

\f{z-Ks)\ < C(1 -I- |z-is|)^e-«|i'"(^"'")l < C(l + s/T)^e^'*/^ V \z\ = t. 

Since B/t < 1, it follows that |/(z~^s)| is absolutely integrable against e^'^ds. The 
Fubini-Tonelli theorem then tells us that 

/(^)dP(^) = J_I (i) 9tl^dz, (39) 



2TTi J\z\=T \^ J ^ 

where 

poo 

<!>j:{z) = / e-'f{zs)ds. (40) 

Here, we will not be concerned with the numerical evaluation of ^f{z). This depends 
upon the choice of / and, in most cases, the integral can be performed explicitly. As we 
have already mentioned in the second section, for orthogonal polynomials, the Gauss- 
Laguerre quadrature rule is more appropriate. The function ^f{z) is analytic in the 
strip defined by the equation |Im(z)| < 1/B. The assertion follows from Proposition ^ 
by noticing that ^f{z) is a function of the type defined by Eq. (|26|l . with the measure 
dP{uj) obtained from the density e~'^, if w > 0, and 0, otherwise. For this measure, 
To = 1, whereas the correlation (characteristic) function is given by the well known 
expression C{t) = 1/(1 — it). Eq. then gives an alternative way to evaluate ^f{z), 
namely 

'^A^) = r T^dK. (41) 



B 



1 — izK 



We now go back to Eq. H39|l and replace the right-hand side integral with a trapezoidal 
quadrature sum. We obtain 

r 1 ^ 



In the notation of Fornaro, if we let g(z) = f{l/z)C{—iz)/z, then the difference between 
the left-hand side and the right-hand side of Eq. (|42|l is denoted by Epf{g). Notice that 
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the function f{l/z) is analytic in the domain \z\ > B, whereas the function C{—iz) 
is analytic in the domain \z\ < tq. It follows that g{z) is analytic on the annulus 
B < \z\ < Tq. If < p-i < 1 < p2 are chosen such that B < piT < r < p2T < tq, then 
Corollary 5.1 of Ref. |23 asserts that 

If Mi,M2 > are bounds of g(z) on the circles of equations \z\ = piT and \z\ = p2T, 
respectively, then Eq. iBl^ implies 

\EN{g)\ < , + , (44) 

from which the exponential decay to zero follows quickly if one only remembers that 
pi < 1 and p2 > 1. 

The result of Fornaro generalizes previous error estimates of Lyness and Delves^ It 
can produce accurate estimates of the error, provided that the upper bound expressed 
by Eq. H44|l is more carefully established. However, such error estimates are not within 
the scope of the present paper. 



VI. CONCLUSIONS 

We have established that generalized moments of spectral distributions against Paley- 
Wiener type functions can be computed from short-time information via Eqs. I0J and 
However, it is required that the correlation function be defined by continuity on the 
entire complex place with the possible exception of two branch cuts on the imaginary 
axis. The values outside the small open disk about the origin do not have to be known, 
yet their existence is important. Our effort to prove Theorem |21 for distributions was 
motivated by the fact that we wanted to maintain the polynomials and the complex 
exponentials of the form exp(izt) in the class of admissible functions. Indeed, these 
functions are of Paley- Wiener type because they are the Fourier-Laplace transform of 
Dirac's delta function or its distributional derivatives. 

In practical applications, we will probably be more interested in the case when the 
distributions are represented by some integrable functions /(k) vanishing outside some 
bounded interval [—B, B]. If the function is also infinitely differentiable on [—B, B], say, 
if it is 

f(K) = I ''''P [-'^'/(S" - K^)] , if |k| < B, 
^ ' \ 0, otherwise, ^ ^ 

the original version of the Paley- Wiener theorem tells us that, for each N = 1,2, . . ., 
there is a constant Cn > such that 

l/WI <C^(l + k|)~'^e^|i-(^)l. 

This says that f{z) decays faster than any polynomial on the real axis. Although, they 
cannot compete with a Gaussian in terms of decay of the tail, the Paley- Wiener func- 
tions obtained by transforming arbitrarily smooth and compactly supported functions 
are numerically more useful than polynomials or complex exponentials. Due to their 
integrability, they can form delta sequences and be utilized to generate representations 
of the type given by Eq. ®. In this respect, it is useful to remember that, if <t > and 
Wo G K., then /((z — LUo)/a) is a Paley- Wiener function as well. This follows easily from 
the equality 



Again, we mention that if the correlation function cannot be extended analytically 
beyond the strip 5, the integral transformation is still valid provided that the support 
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of the distribution /(k) is included in the interval (— 2to,2to). If the support of the 
distribution /(k) is further restricted to be included in the interval {—tq^tq), then the 
evaluation of the generalized moments f{uj)dP(uj) can be performed exponentially fast 
by means of Eq. (|42() . Although we lose the ability of forming delta sequences, we still 
have the possibility of using Paley- Wiener functions that decay to infinity. Such functions 
are more useful than polynomials in the context of maximum entropy techniques^sSiii 
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